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We investigate the temperature dependence of the Landau gauge gluon and ghost propagators in 
lattice QCD with two flavors of maximally twisted mass fermions. For these propagators we provide 
and analyze data which corresponds to pion mass values between 300 and 500 MeV. For the gluon 
propagator we find that both the longitudinal and transversal component change smoothly in the 
crossover region, while the ghost propagator exhibits only a very weak temperature dependence. 
For momenta between 0.4 and 3.0 GeV we give a parametrization for our lattice data. It may serve 
as input to studies which employ continuum functional methods. 
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I. INTRODUCTION 

Hadronic matter may experience a phase transition or 
crossover from a confined phase, with broken chiral sym- 
metry, to a deconfined chirally symmetric phase. The 
latter phase is characterized by a state called quark gluon 
plasma. It is widely believed this state of matter has been 
passed in the early universe, and efforts are undertaken 
to reproduce it experimentally in heavy-ion collisions as 
in present collider experiments at RHIC (BNL) or with 
the ALICE and CMS detectors at LHC (CERN). 

The possible existence of such a phase transition was 
first discussed [Ij in the context of Hagedorn's thermody- 
namic model [2, y| as a way to evade the consequence of a 
maximal hadronic temperature. First numerical evidence 
has been found already in the early days of lattice QCD 
(LQCD) [4j, 5j. In fact, the lattice formulation offers an 
ab-initio approach to study such aspects of QCD nonper- 
turbatively. This remains true as long as the chemical 
potential is small compared to the temperature. Over 
the last decade immense computational resources were 
dedicated to reach at a consistent picture of QCD at fi- 
nite temperature and zero chemical potential (for recent 
reviews see (l-Q)- 

LQCD, however, is not the only framework to tackle 
nonperturbative problems of QCD at zero or non-zero 
temperature. Powerful continuum functional methods 
exist as well, like for instance in the context of Dyson- 
Schwinger (DS) equations [Tol - [l3| or functional renor- 
malization group (FRG) equations [H [Tsj i. which also 
allow to address such problems. 

Within these frameworks (see, e.g., the recent re- 
views [H, the Landau gluon and ghost propaga- 
tors appear — together with the corresponding vertices — 
as the main building blocks in the formulation of the DS 
or FRG equations; they constitute part of the solutions 
of the latter. These functional methods though come 
with a potential source of error: To solve the (infinite) 



tower of equations it has to be truncated appropriately, 
and so the solutions beyond the far-infrared regime, de- 
pend on how the truncations are done, in particular in 
the momentum range around 0(1) GeV. Therefore, in- 
dependent information, at best from first principles, is 
welcome to improve these (unavoidable) truncations. 

LQCD calculations allow to provide the Landau gluon 
and ghost propagators in an ab-initio way. The available 
momentum range, however, is restricted from above by 
the lattice spacing and from below by the available lattice 
volume (up to a further uncertainty related to so-called 
Gribov copies [H, HI])- Despite these restrictions, an 
impressive amount of data has been produced over the 



last years for these propagators at zero (see, e.g., 20 
21 1 and references therein) and non-zero temperature [19 
20 , [22l - [35l | . This data has allowed for a variety of cross- 



checks with corresponding results from the continuum 
functional methods. 

Conventional lattice calculations of QCD at finite tem- 
perature often employ the Polyakov loop and the chiral 
condensate to probe for the (de)confinement and chiral 
phase transition, respectively. In recent years it turned 
out that these and other observables can be also calcu- 
lated from DS and FRG equations, where the Landau 
gauge gluon and ghost propagators are used as input 
29, 36 42]. Even an extrapolation to the notoriously dif- 
ficult re gim e of nonzero chemical potential seems to be 
possible |43l.l43|. Appropriate lattice data for the propa- 
gators close to the continuum limit is therefore essential 
to assist those efforts. 

In a recent study [32| we have provided such data for 
quenched QCD (see, e.g., [i^ for a first application). We 
could show that the longitudinal (electric) component of 
the gluon propagator may be used to probe the thermal 
phase transition of pure SU{3) gauge theory. Similar was 
shown in [ssj. The present article further complements 
the available lattice data with finite-temperature data 
for the Landau gauge gluon and ghost propagators from 
LQCD with Nf = 2 dynamical fermion flavors. To the 
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best of our knowledge, there are only two studies which 
have provided such data in the past [s^, • 

For our study we adopt a lattice formulation that em- 
ploys the Symanzik-improved gauge action for the glu- 
onic field and the twisted mass Wilson-fermion action 
for the fermionic part. The latter ensures an automatic 
0(a) improvement provided the Wilson K-parameter is 
tuned to maximal twist (for further details we refer to 
p7l lisj). Our data is based on the vast set of gauge 
field configurations that has been generated by the tmfT 
Collaboration. The tmfT Collaboration has^xplored the 
complicated phase structure of the theory i49;j and is still 
investigating the smooth crossover region from the con- 
fining and chirally broken regime at low temperature T 
to the deconfinement and chirally restored phase at high 
T [50I, 'Sl'l . These studies are restricted so far to pseudo- 
scalar meson (pion) masses from 300 MeV to 500 MeV. 
To set the scale we use results at T = of the Euro- 
pean Twisted Mass (ETM) Collaboration Specifi- 
cally, our data for the gluon and ghost propagators cov- 
ers the whole crossover regime at three pion mass val- 
ues: TO^ ~ 316, 398 and 469 MeV. At these values the 
crossover regime is characterized by a very smooth be- 
havior of the chiral condensate and the Polyakov loop as 
well as their susceptibilities. Moreover, one observes for 
these settings the breakdown of chiral symmetry and the 
deconfinement phase transition occur at slightly differ- 
ent temperatures = and Tdccont, respectively, in 
agreement with observations reported in j53l|.^ 

The paper is organized as follows. In Section we 
give all lattice parameters and outline the setup of our 
Monte Carlo simulations. Section IIIII recalls the defini- 
tions of the gluon and ghost propagators on the lattice 
in the Landau gauge. Data and fits for various temper- 
ature and pion mass values are presented in Section IIVI 
Conclusions are drawn in Section IVl 



II. LATTICE ACTION AND SIMULATION 
PARAMETERS 



Our study is based on gauge field configurations pro- 
vided by the tmfT Collaboration. These configurations 
were generated on a four-dimensional periodic lattice of 
spatial linear size N^r ~ 32 and a temporal extent of 
Nr = 12 for a mass-degenerate doublet of twisted mass 
fermions, cf. the review in Ref. j47| . The corresponding 
gauge action is the tree-level Symanzik improved action 



^ Detailed results are presented in a recent update to |51| . We 
thank the tmfT Collaboration for providing us their Tc-data prior 
to publication. 
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with f3 — 2Nc/gQ , ci = — , cq = 1 — 8 ci and go being 
the bare coupling constant. f/^^J represents quadratic 
and U]^^^ rectangular Wilson loops built from the link 
variables Uxfj. € SU{3). One important feature of this 
gauge action is its inherent 0{a) improvement. For more 
details see Refs. [13, HH]- The Wilson fermion action with 
an additional parity-flavor symmetry violating improve- 
ment term reads [56| 

S'_f[J7, '0, V'] = 

^X{x) (1 - kDw[U] + 2iKa^o75T3) x(x) , (3) 

where the Pauli matrix T3 acts in flavor space and the 
fermionic fields are expressed in terms of the twisted basis 
which is related to the physical fields by 

■0 = -^(1 + 175x3)^ and V = X-^(1 + 1757-3) . (4) 
The Wilson covariant derivative acts on these as 
DwPWx) = ^((r-7^)[/^(x)0(x-f A) 

{r + j,)Ul{x~fi)^P{x-fLj), (5) 

and the quark mass is set by the twisted mass parameter 
fiQ and the hopping parameter k = (2amQ + 8r)^^, pa- 
rameterizing the untwisted bare quark mass component. 
Here a is the lattice spacing and r — 1. Note that for 
any finite /3 the value for k gets corrections through mass 
renormalization. 

It must be noted that maximal twist is accomplished 
by tuning the hopping parameter to its critical value Kc, 
where the untwisted theory would become massless. At 
maximal twist one achieves an automatic 0(a)-improved 
fermion formulation (STj . The hopping parameter enter- 
ing the simulation is based on a set of /3- values for which 
KciP) was provided by the ETM Collaboration [sll- For 
intermediate l3 the corresponding Kc values are obtained 
through an interpolation as described in Ref. [51.]. The 
bare twisted mass parameter a/io has been adjusted such 
as to keep the physical pion mass constant along our 
scans in the bare inverse coupling /3. 

As usual in finite temperature QCD, the imaginary 
time extent Nr corresponds to the inverse temperature 
T^^ — Nra{f3). To quote it in physical units we use 
interpolated — as well as slightly extrapolated — data for 
the lattice spacing reported for (3 = 3.90, 4.05 and 4.20 
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by the ETM Collaboration [13 (see also Ref. [U). We 
restrict our analysis to lattice spacings a < 0.09 fm. 

For the reader's convenience all parameters, like /3, 
the corresponding lattice spacings, temperatures, pion 
masses, the number of independent configurations and 
other relevant values are collected in TableU] Note that 
there, for definiteness, "pion mass" corresponds to the 
charged pion. In Tablellll we provide also the respective 
pseudo-critical couplings (3c and the corresponding tem- 
peratures and Tdoconf for the three pion mass values we 
use. These temperatures were obtained from fits around 
the maxima of the chiral susceptibility and from the 

behavior of the (renormalized) Polyakov loop (*He(L))^, 
respectively (see the revised version of Ref. [51]). 

III. THE GLUON AND GHOST 
PROPAGATORS ON THE LATTICE 

Gluon and ghost propagators are gauge dependent 
quantities. As in Ref. [32] we focus on Landau gauge 
and therefore have to transform the (unfixed) tmfT gauge 
ensemble until it satisfies the corresponding gauge con- 
dition. In differential form it reads 

4 

V^A^ - J2 iM^ + A/2) - A^ix - fi/2)) = (6) 
with the lattice gauge potentials 

Af,{x + fl/2) ^ —^{U^f, - ) \traceless ■ (7) 

2iago ^ 

To render the link variables satisfying this condition one 
maximizes the gauge functional 

Fu[g]^^Y.'^cTT{glU,,g,+,,) (8) 

by successive local gauge transformations gx acting on 
the link variables as follows 

f/.M ^ = 9lUx^9.+^. , 9. e SU{3) . (9) 

In order to achieve this, we subsequently apply two 
methods, first simulated annealing (SA) and then over- 
relaxation (OR). OR is applied to finally satisfy the 
gauge condition (Eq. ([6])) with a local accuracy of 

maxfHeTr[V^A,^V,Ay < IQ-^^ ^ (10) 

while SA to reduce the Gribov ambiguity of lattice Lan- 
dau gauge, by favoring gauge-fixed (Gribov) copies with 
large values for -F'[/[5], see Jfl, 21, 59-64]. For this the 
SA algorithm generates gauge transformations {gx} ran- 
domly by a Monte Carlo chain with a statistical weight 
^ exp{Fij[g]/Tsa)- The "simulated annealing tempera- 
ture" Tsa is a technical parameter which is monotonously 
lowered. Our annealing schedule is specified by a hot 



start at T^a = 0.45, after which Tga is continuously low- 
ered in equal steps until Tga = 0.01 is reached. We ap- 
ply 3500 SA steps between these two temperatures, and, 
for better performance, also added a few microcanonical 
steps to each (heatbath) step. Since we apply large num- 
ber of SA steps to maximize Fjj [g] , we restrict ourselves 
to one Gribov copy per configuraton. 

Our first quantity of interest is the gluon propagator 
defined in momentum space as the ensemble average 

DfM = (A;ik)Ali-k)) , (11) 

where (• • • ) represents the average over configurations. 
A^(fc) denotes the Fourier transform of the gauge poten- 
tial ([7]) and fc^t € {—Nf^/2, N^I2\ is the lattice momentum 
(/i = 1, . . . , 4), which relates to the physical momentum 
as 

Henceforth we will use the notation [N^; Nr) = {Ni \ N4) 
where i = 1,2, 3. 

For non-zero temperature Euclidean invariance is bro- 
ken, and it is useful to split D'^^{q) into two components, 
the transversal Dt ( "chromomagnetic" ) and the longitu- 
dinal Dl ( "chromoelectric" ) propagator, respectively, 

DfM^6'^'{P^Mqlq^) + P;:,DL{ql,q^)) . (13) 

The fourth momentum component q^ conjugate to the 
Euclidean time (Matsubara frequency) will be restricted 
to zero later on. For Landau gauge Pjj^ represent pro- 
jectors transversal and longitudinal relative to the time- 
direction (/i = 4): 

P^. = (1 - 5,,){1 - 5,,) (^S,, ~ , (14) 

For the propagators Dt.l [or their respective dimension- 
less dressing functions Zx.Lil) — 'Z^-Dt,l('z)] we find 

Driq) = ^ ^^I^(fc)4^(-fc) - ^At{k)Al{-k)j 

(16) 

and 

DUq) = ]^ (1 + |l) {AUk)Ati-k)) , (17) 

where Ng — — 1 and Nc — 3. The zero-momentum 
propagator values are then defined as 

^t{0) - 5^E(^°(0M?(0)), (18) 

^ i— 1 

DUO) = ^(I^(O)I^(O)) . (19) 
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m„ [MeV] 


T [McV] 


a [fm] 


raja 


ro ■ T 




316(16) 


187 


8.77(47) 


10" 


2 


4.81 


0.40 


293 


316(16) 


199 


8.25(22) 


10" 


2 


5.17 


0.43 


299 


316(16) 


215 


7.65(13) 


10" 


2 


5.63 


0.47 


255 


316(16) 


222 


7.39(12) 


10" 




5.84 


0.49 


273 


316(16) 


225 


7.31(12) 


10" 




5.91 


0.49 


151 


316(16) 


228 


7.22(11) 


10" 




5.98 


0.50 


250 


316(16) 


230 


7.14(11) 


10" 




6.05 


0.50 


113 


316(16) 


235 


6.98(11) 


10" 




6.19 


0.52 


290 


398(16) 


193 


8.51(32) 


10" 


2 


4.99 


0.42 


159 


398(16) 


199 


8.25(22) 


10" 


2 


5.17 


0.43 


173 


398(16) 


215 


7.65(13) 


10" 


2 


5.63 


0.47 


209 


398(16) 


228 


7.20(11) 


10" 


2 


6.00 


0.50 


198 


398(20) 


236 


6.98(11) 


10" 




6.19 


0.52 


156 


398(16) 


241 


6.82(10) 


10" 




6.34 


0.53 


150 


398(20) 


246 


6.69(10) 


10" 




6.46 


0.54 


271 


398(20) 


248 


6.62(10) 


10" 




6.53 


0.54 


226 


398(20) 


254 


6.47(10) 


10" 




6.68 


0.56 


113 


469(24) 


222 


7.42(12) 


10" 


2 


5.81 


0.48 


146 


469(24) 


228 


7.20(11) 


10" 


2 


6.00 


0.50 


348 


469(24) 


235 


6.98(11) 


10" 


2 


6.19 


0.52 


120 


469(24) 


243 


6.77(10) 


10" 


2 


6.39 


0.53 


210 


469(24) 


247 


6.67(10) 


10" 


2 


6.48 


0.54 


250 


469(24) 


250 


6.57(10) 


10" 


2 


6.58 


0.55 


256 


469(24) 


254 


6.47(10) 


10" 


2 


6.68 


0.56 


152 


469(24) 


258 


6.38(10) 


10" 


2 


6.78 


0.56 


150 


469(24) 


266 


6.19(12) 


10" 


2 


6.98 


0.58 


200 





K,(/3) 


a ■ fj-oiP) 


3.8400 


0.162731 


0.00391 


3.8800 


0.161457 


0.00360 


3.9300 


0.159998 


0.00346 


3.9525 


0.159385 


0.00335 


3.9600 


0.159187 


0.00331 


3.9675 


0.158991 


0.00328 


3.9750 


0.158798 


0.00325 


3.9900 


0.158421 


0.00319 


3.8600 


0.162081 


0.00617 


3.8800 


0.161457 


0.00600 


3.9300 


0.159998 


0.00561 


3.9700 


0.158927 


0.00531 


3.9900 


0.158421 


0.00517 


4.0050 


0.158053 


0.00506 


4.0175 


0.157755 


0.00498 


4.0250 


0.157579 


0.00493 


4.0400 


0.157235 


0.00483 


3.9500 


0.159452 


0.00779 


3.9700 


0.158926 


0.00752 


3.9900 


0.158421 


0.00738 


4.0100 


0.157933 


0.00718 


4.0200 


0.157696 


0.00708 


4.0300 


0.157463 


0.00699 


4.0400 


0.157235 


0.00689 


4.0500 


0.157010 


0.00680 


4.0700 


0.156573 


0.00662 



Zt 


Zl 


Zj 


0.6380(80) 


0.6264(108) 


0.66862(32) 


0.6208(34) 


0.6139(50) 


0.66939(10) 


0.6117(70) 


0.6116(103) 


0.67264(14) 


0.6156(59) 


0.6122(84) 


0.67252(15) 


0.6151(71) 


0.6101(102) 




0.6120(63) 


0.6079(114) 


0.67388(14) 


0.6146(77) 


0.5982(126) 




0.6092(71) 


0.6118(97) 


0.67536(14) 


0.6191(64) 


0.6147(98) 


0.66730(21) 


0.6202(54) 


0.6192(76) 




0.6076(59) 


0.6080(84) 


0.67005(23) 


0.6087(59) 


0.6119(89) 




0.6075(139) 


0.6090(231) 


0.67293(13) 


0.6156(87) 


0.6112(115) 


0.6784(60) 


0.6036(78) 


0.6063(101) 


0.67504(16) 


0.6001(64) 


0.6005(94) 


0.67517(11) 


0.6029(119) 


0.6167(178) 


0.67519(26) 


0.6121(55) 


0.6020(82) 


0.67142(16) 


0.6116(80) 


0.6024(110) 


0.67128(14) 


0.6098(70) 


0.6041(102) 


0.67271(22) 


0.6086(54) 


0.6093(73) 


0.67322(12) 


0.5947(54) 


0.5927(77) 


0.67147(14) 


0.6013(72) 


0.6017(101) 


0.67388(14) 


0.6033(80) 


0.6028(121) 


0.67353(16) 


0.5971(62) 


0.6072(91) 


0.67485(17) 


0.5972(143) 


0.6119(195) 


0.67829(32) 



TABLE I: The pion mass values, the temperature T, both in MeV, the lattice spacing a in fm, the chirally extrapolated Sommer 
scale ro [1^, ro T, and the number riconf of independent configurations used for the analysis are shown in the left subtable for 
all simulation ensembles. 

In the central subtable the values of the inverse bare coupling /? used in the simulations, the critical hopping parameter Kc and 
the bare twisted mass are additionally shown. The spatial {Ncr ~ 32) and temporal {Nt — 12) sizes are the same for all 
ensembles. The number ricopy of gauge copies was fixed to 1. 

In the right subtable the renormalization factors for the transverse and longitudinal gluon as well as for the ghost dressing 
function denoted by Zt, Zl and Zj, respectively (see the text), are given for the renormalization scale /i — 2.5 GeV. 



tmfr ensemble 


A12 


B12 


C12 


m„ [MeV] 


316(16) 


398(20) 


469(24) 


/3c from 


3.89(3) 


3.93(2) 


3.97(3) 


[MeV] 


202(7) 


217(5) 


229(5) 


/3c from (5Re(L))^ 




4.027(14) 


4.050(15) 


Tdcconf [MeV] 




249(5) 


258(5) 



TABLE II: Extracted (pseudo-) critical couplings Pc and the 
crossover-temperatures and Tdcconf for the three ensembles 
A12, B12, and C12. Corresponding pion masses (from [5l|]) 
are listed in the first row. Ensemble names indicate Nr = 12. 



Note that we have neglected a possible 0{a) improvement 
related to the use of the improved gauge action Eq. ^ . 
The Landau gauge ghost propagator is given by 



2-ni(k/N)-(x-y) [j\//-l]ob^ 



5"' J{q)lq' 



(20) 



where the four- vector (k/N) = {k^/N^). J{q) denotes 
the ghost dressing function. The matrix M is the lattice 



Faddeev-Popov operator 



xy 



El 



(21) 



with 



<^ = $HeTr[{r^T''}([/,,^ + C/,_^,^)] 
B2;'^ = 2-5HeTr[r''T'^C/,J, 
Cf = 2 • mcTT[T''T''a 



Ng, i.e. the Hermitian 



written in terms of T", a = 1, 
generators of the su(3) Lie algebra normalized accord- 
ing to Tr[rT''] = (5''V2- For the inversion of M we 
use the pre-conditioned conjugate gradient algorithm of 
[B^ I with plane-wave sources ■0c with color and position 
components = 6xp(27rifc • (x/N)). 

To reduce lattice artifacts, we apply cylinder and cone 
cuts to our data [6^1 • Specifically we consider only di- 
agonal and slightly off-diagonal momenta for the gluon 
propagator and diagonal momenta for the ghost propa- 
gator. Moreover, only modes with zero Matsubara fre- 
quency (fc4 = 0) are used. 
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IV. GLUON AND GHOST PROPAGATOR 
RESULTS 

A. Momentum dependence 

Data for the unrcnormalized transverse (Zt) and lon- 
gitudinal (Zl) gluon dressing functions and also for the 
ghost dressing function (J) is shown in Fig.[TJ We show 
it versus the physical momentum q = \q\ for selected 
temperatures and for three pion masses (panels from 
top to bottom are for ~ 316, 398 and 469 MeV, re- 
spectively). The corresponding renormalized functions, 
in momentum subtraction (MOM) schemes, can be ob- 
tained from 

Z^j.2{q,^l) EE Zt,l(m)^T,l(<z), 

J-"(g,A*) EE Zj{lJi)J{q) (22) 

with the Z- factors being defined such that ZJj?'^l{^, ji) — 
J^'^^{^,fi) = 1. For a renormalization scale oi fi — 
2.5 GeV the Z-factors are quoted in TableU] 

Fig. m also shows curves connecting data points of same 
temperature. These were obtained from fits to the data 
for momenta 0.4 GeV < q < 3.0 GeV. These curves may 
serve as input to studies of the corresponding DS or FRG 
equations. 

More specifically, for the gluon dressing function we 
employed (analogously to our quenched study ^) the 
Gribov-Stingl formula ^ 113 

2 c(l + rfg^") 
^fit(<z) - q (^2^.^2)2 ^^2 ' (23) 

which has been also used in [3l|, and appears in 
the context of the so-called "Refined Gribov-Zwanziger" 
framework [g^, [tOI ■ But we found it sufficient to set = 
and n = 1. This fits well the data for q e [0.4, 3.0] GeV 
and gives excellent Xdof values. The latter together with 
the results for the fit parameters are listed in Table lllll 
Note again that Fig.[T] shows data and the correspond- 
ing fits only for a selected range of temperatures, but 
Table Um gives the fit parameters for all available tem- 
peratures. We cannot exclude that the 5^-term is needed 
for smaller momenta. If true, it would indicate the oc- 
currence of a pair of complex-conjugate poles. 

For momenta above 3 GeV the fit fails to describe the 
data. In this range logarithmic corrections are expected 
to become important. 

For the ghost dressing function we propose to use a fit 
formula like 

Jfit{q)= (^^) +h (24) 

In a first attempt we also tried hq'^/ {q^ + m'^h) for the last 
term with mgh as a free parameter, but this was always 
found being consistent with mgh — 0. We therefore omit 
such infrared mass parameter and only keep a constant 
term h in the ultraviolet limit. 



Fit results for the fitting range [0.4 GeV, 4.0 GeV] are 
presented also in Table lllll One notes that our Xdof values 
are far from being optimal, in particular for the lower 
temperatures. Deviations typically occur at the lowest 
momenta. But this could not be cured, e.g., by a mass 
term m^^ alone. However note, the maximal deviations 
of fit and data points do not exceed 5%. 



B. Temperature dependence 

We now look at the temperature dependence of the 
dressing functions, where our temperature values cover 
the chiral restoration and the deconfinement phase tran- 
sition, with the latter being signaled by a peak in the 
Polyakov loop susceptibility. These two crossover phe- 
nomena typically occur at different temperatures and will 
be denoted and Tdeconf, respectively, in what follows 
(see Tableini for their values). 

Looking once again at Fig.[Tl we see that the momen- 
tum dependences of Z^lq) and Zxiq) change differently 
with temperature, irrespective of m,r- In fact, while the 
(unrenormalized) transverse dressing function Zt seems 
to be relatively insensitive to the temperature, the curves 
describing Z^i^q) fan out for momenta below the renor- 
malization scale /i = 2.5 GeV. A stronger temperature 
dependence we also observed for DL{q) for pure SU{i) 
gauge theory [s^ ]. though there it was found to be much 
more pronounced due to the existence of a first order 
phase transition [s^ . 

These observations are seen more clearly in Fig. [21 
where we show ratios of the renormalized dressing func- 
tions or propagators 

RrAq^T) = D^T%q,T)/D^T%q,%^^), (25) 
RG{q,T) = G"^™(g,r)/G^^"(g,r,„i„) (26) 

as functions of the temperature T for 6 fixed (interpo- 
lated) momentum values g 7^ 0, and for different pion 
masses (panels from top to bottom). For better visibil- 
ity, ratios are normalized with respect to the respective 
left-most shown temperature in Fig. [51 

Looking at Fig.[21 we see Rhiq, T) decrease more or less 
monotonously with temperature in the crossover region, 
and this decrease is stronger the smaller the momentum. 
Rxiq, T) instead signals a slight increase within the same 
range, and the ghost propagator (at fixed low momenta) 
seems to rise a bit around T ~ Tdeconf- 

Fig. [21 does not show ratios at zero momentum, but 
we show it for Rt in Fig. [31 (upper row), again versus 
T and from left to right for different pion masses. For 
TOtt ~ 398 MeV, for example, Rt{0,T) clearly rises to- 
wards Tdeconf, whereas there are only weak indications 
for such a behavior for Rt{0,T) for the other two data 
sets. Much more statistics is necessary to resolve that. 

The lower panels of Fig. [31 show data for the inverse 
renormalized longitudinal propagator D£'^" at zero mo- 
mentum, again versus temperature and from left to right 
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Zt fits 


13 




djo? 


a r 


xlaf 


3.8400 


1.334(132) 


0.744(138) 


0.415(14) 


0.19 


3.8800 


1.183(59) 


0.872(79) 


0.390(7) 


0.16 


3.9300 


1.032(101) 


1.013(188) 


0.370(13) 


0.12 


3.9525 


1.049(93) 


0.941(159) 


0.370(11) 


0.16 


3.9600 


0.932(103) 


1.148(233) 


0.355(14) 


0.13 


3.9675 


1.023(108) 


0.979(206) 


0.369(13) 


0.27 


3.9750 


0.938(119) 


1.165(277) 


0.359(15) 


0.09 


3.9900 


0.914(98) 


1.143(224) 


0.358(13) 


0.38 


3.8600 


1.271(116) 


0.799(143) 


0.401(13) 


0.19 


3.8800 


1.218(83) 


0.808(107) 


0.391(9) 


0.11 


3.9300 


0.982(91) 


1.092(177) 


0.356(12) 


0.11 


3.9700 


1.034(86) 


0.928(160) 


0.366(10) 


0.67 


3.9900 


1.171(265) 


0.699(430) 


0.383(42) 


0.21 


4.0050 


1.006(145) 


0.924(281) 


0.366(19) 


0.15 


4.0175 


0.922(133) 


1.087(297) 


0.356(18) 


0.12 


4.0250 


0.845(94) 


1.278(270) 


0.346(13) 


0.11 


4.0400 


0.824(173) 


1.252(502) 


0.348(24) 


0.08 


3.9500 


0.952(78) 


1.127(169) 


0.348(10) 


0.35 


3.9700 


0.993(100) 


1.003(199) 


0.351(12) 


0.32 


3.9900 


0.883(96) 


1.203(236) 


0.343(13) 


0.19 


4.0100 


0.998(82) 


0.945(157) 


0.368(10) 


0.23 


4.0200 


1.010(80) 


0.918(156) 


0.356(9) 


1.08 


4.0300 


0.900(87) 


1.108(223) 


0.349(11) 


0.67 


4.0400 


0.870(103) 


1.131(275) 


0.340(14) 


0.25 


4.0500 


0.826(89) 


1.251(260) 


0.344(13) 


0.25 


4.0700 


0.971(250) 


0.894(559) 


0.371(32) 


0.01 





Zl fit 


s 




c/a^ 


d/a^ 


a r 


xLf 


1.868(142) 


0.420(78) 


0.510(12) 


0.13 


1.729(60) 


0.463(40) 


0.486(5) 


0.66 


1.371(104) 


0.647(111) 


0.431(10) 


0.36 


1.218(76) 


0.757(96) 


0.411(8) 


0.22 


1.208(88) 


0.744(122) 


0.405(10) 


0.05 


1.256(83) 


0.682(107) 


0.410(9) 


0.10 


1.099(93) 


0.852(156) 


0.387(11) 


0.07 


0.982(70) 


1.026(143) 


0.369(9) 


0.38 


1.811(109) 


0.459(66) 


0.501(9) 


0.32 


1.637(88) 


0.523(63) 


0.477(8) 


0.66 


1.389(84) 


0.651(86) 


0.437(8) 


0.37 


1.151(71) 


0.819(110) 


0.397(8) 


0.48 


1.035(171) 


0.965(310) 


0.383(21) 


0.25 


1.092(104) 


0.784(162) 


0.381(12) 


0.25 


1.035(87) 


0.891(164) 


0.373(10) 


0.44 


1.020(78) 


0.904(144) 


0.368(9) 


0.21 


0.837(111) 


1.237(329) 


0.340(15) 


0.12 


1.338(80) 


0.645(89) 


0.427(8) 


0.63 


1.290(104) 


0.648(127) 


0.416(11) 


0.14 


1.126(91) 


0.797(139) 


0.389(10) 


0.28 


1.009(61) 


0.937(115) 


0.371(7) 


0.32 


1.027(67) 


0.946(124) 


0.375(8) 


0.07 


1.096(79) 


0.775(135) 


0.380(9) 


0.48 


0.941(84) 


1.018(191) 


0.359(11) 


0.22 


1.053(69) 


0.781(126) 


0.367(8) 


0.67 


0.792(130) 


1.291(405) 


0.328(34) 


0.15 



Zj fits 






k 


2 

Xdaf 


0.41822(7) 
0.37046(9) 
0.35672(9) 


1 nQ16('61'l 

1.0904(19) 
1.1355(39) 
1.1387(36) 


5111 ('781 

0.4950(23) 
0.5438(55) 
0.5462(52) 


69 
8.59 
2.29 
2.40 


0.34636(7) 


1.1501(33) 


0.5642(54) 


5.20 


0.33093(8) 


1.1571(30) 


0.5736(56) 


7.01 


0.4464(21) 


1.0419(40) 


0.4444(35) 


31.1 


0.3802(19) 


1.0962(53) 


0.4945(58) 


21.4 


0.3380(07) 
0.441(84) 
0.3189(09) 
0.3155(08) 
0.3135(21) 


1.1466(25) 
0.90(14) 
1.1623(31) 
1.1630(25) 
1.2200(90) 


0.5612(39) 

0.41(06) 
0.5834(60) 
0.5853(52) 
0.696(21) 


24.0 
0.21 
6.4 
17.4 
0.93 


0.3621(09) 
0.3519(07) 
0.3405(12) 
0.3295(07) 
0.3270(08) 
0.3130(08) 
0.3100(11) 
0.3082(11) 
0.2913(20) 


1.1325(38) 
1.1426(31) 
1.1555(45) 
1.1607(24) 
1.1557(27) 
1.1538(26) 
1.2015(55) 
1.2109(51) 
1.2156(79) 


0.5397(53) 
0.5554(48) 
0.5787(81) 
0.5873(48) 
0.5778(51) 
0.5660(48) 
0.648(12) 
0.677(12) 
0.682(19) 


2.3 
4.2 
2.3 
13.1 
7.0 
7.9 
2.3 
0.94 
0.39 



TABLE III: Results from fits with the Gribov-Stingl formula Eq. (|23|) for the unrenormalized Zt (left subtable) and Zl (center 
subtable) gluon dressing functions. The fit range is [0.4 : 3.0] GeV. The values in parentheses indicate the fit errors estimated 
with the bootstrap method. The parameters h and n were fixed to 6 = and n = 1, respectively. 

Fit results for the unrenormalized ghost dressing function Zj with the fitting function according to Eq. (|24p are presented in the 
right subtable. The momentum fitting ranges are here [0.4 : 4.0] GeV. The pion mass values are = 316(16) MeV (upper), 
= 398(20) MeV (middle) and = 469(24) MeV (bottom subtables), respectively. 



for different pion masses. This quantity can be identified 
with an infrared gluon screening mass, and we clearly see 
it to rise with temperature in the crossover region. This 
again shows that this infrared gluon screening mass may 
serve as an useful indicator for the finite-temperature 
crossover of the quark-gluon system. However, we should 
keep in mind that the zero-momentum results for the 
gluon propagators are influenced also by strong finite- 
size and Gribov copy effects, which we could not analyze 
here. 



V. CONCLUSION 

We have presented data for the Landau gauge gluon 
and ghost propagators for lattice QCD at finite temper- 
ature with Nf = 2 twisted mass fermion fiavors. Our 
data is for a momentum range of 0.4 GeV to 4.0 GeV 
and was obtained on gauge field configurations produced 
by the tmfl Collaboration. This has allowed us to ex- 
plore the propagator's momentum dependence over the 
whole temperature range of the crossover region, and this 
separately for three (charged) pion mass values between 



300 MeV and 500 MeV. We find that the propagators 
change smoothly passing through the crossover region 
and the most significant change is seen for the longitudi- 
nal (i.e., electric) component of the gluon propagator. 



We also provide fitting functions for our data. These 
and the corresponding fit parameters, given in Table lllll 
may serve as interpolation functions of our data when 
used as input to studies which employ continuum func- 
tional methods to address problems of QCD at finite tem- 
perature. Actually, for both the transversal and longitu- 
dinal gluon dressing function these interpolation func- 
tions give quite a good description of our data for all 
temperatures. For the ghost dressing function, this is 
achieved only for selected temperatures (see Table lllll for 
details). 



We hope our results will help these (continuum-based) 
studies to get further predictions of the behavior of 
hadronic matter close to the transition region that would 
be too difficult if addressed on the lattice directly. 
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FIG. 1: The unrenormalized dressing functions for the transverse gluon Zt (left panel), for the longitudinal gluon Zl (middle 
panel) and for the ghost dressing function J (right panel) are shown as functions of the momentum q [GeV] for different (inverse) 
coupling values /? (i.e. different temperatures) given in the legend. The corresponding pion mass values (from top to bottom 
panels) are m-^ — 316, 398, and 469 MeV. 
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RT{q,T) Rl{clT) Rciq.T) 




T [MeV] T [MeV] T [MeV] 

FIG. 2: Ratios Rt,Rl and Ra for the renormalized transverse DJ?" (left panel), longitudinal Z)£'^" (middle panel) and ghost 
Qren ('j.jgjj^ panel) propagators, respectively, as functions of the temperature T at a few non-zero momentum values q (indicated 
in units of [GeV]. The corresponding pion masses (from top to bottom) are ~ 316, 398 and 469 MeV. The vertical bands 
indicate the chiral and deconfinement pseudo-critical temperatures with their uncertainties (see Table[II)l . 
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FIG. 3: The upper row shows the ratio Rt at zero momentum for the three pion mass values indicated. The lower panels show 
the inverse renormalized longitudinal gluon propagator (D£''")~^ at zero momentum. 



